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ABSTRACT 



1. Introduction 



We review theoretical ideas, problems and implications of neutrino masses and 
mixing angles. We give a general discussion of schemes with three light neutrinos. 
Several specific examples arc analyzed in some detail, particularly those that can 
be embedded into grand unified theories. 



There is by now convinc ing evidence, from the experimental study of atmospheric 
and solar neutrinos l -*- ! ^ ! 4 ^ for the existence of at least two distinct frequencies of 
neutrino oscillations. This in turn implies non- vanishing neutrino masses and a mixing 
matrix, in analogy with the quark sector and the CKM matrix. So apriori the study 
of masses and mixings in the lepton sector should be considered at least as important 
as that in the quark sector. But actually there are a number of features that make 
neutrinos especially interesting. In fact the smallness of neutrino masses is probably 
related to the fact that v's are completely neutral (i.e. they carry no charge which is 
exactly conserved) and are Majorana particles with masses inversely proportional to 
the large scale where lepton number (L) conservation is violated. Majorana masses 
can arise from the see-saw mechanism^*, in which case there is some relation with 
the Dirac masses, or from higher-dimensional non-renormalizable operators which 
come from a different sector of the lagrangian density than any other fermion mass 
terms. The relation with L non-conservation and the fact that the observed neutrino 
oscillation frequencies are well compatible with a large scale for L non-conservation, 
points to a tantalizing connection with Grand Unified Theories (GUT's). So neutrino 



masses and mixings can represent a probe into the physics at GUT energy scales and 
offer a different perspective on the problem of flavour and the origin of fermion masses. 
There are also direct connections with important issues in astrophysics and cosmology 
as for example baryogenesis through leptogenesis ® and the possibly non-negligible 
contribution of neutrinos to hot dark matter in the Universe. 

Recently there have been new important experimental results that have consid- 
erably improved our knowledge. The SNO experiment has confirmed that the solar 
neutrino deficit is due to neutrino oscillations and not to a flaw in our modeling of the 
sun El : the total neutrino flux is in agreement with the solar model but only about 
one third arrives on Earth as v e while the remaining part consists of other kinds of 
active neutrinos, presumably and v T . The allowed amount of sterile neutrinos is 
strongly constrained. The KamLAND experiment has established that T71 from reac- 
tors show oscillations over an average distance of about 180 Km which are perfectly 
compatible with the frequency and mixing angle corresponding to one of the solutions 
of the solar neutrino problem (the Large Angle (LA) solution) Thus the results 
from solar neutrinos have been reproduced and improved by a terrestrial experiment. 
Also the coincidence of the frequency for neutrinos from the sun and for antineutrinos 
from reactors is consistent with the validity of CPT invariance. The validity of this 
symmetry had been questioned because of the puzzling LSND claim of a signal that 
could indicate a third distinct oscillation frequency (hence implying either more than 
three light neutrinos or CPT violation). In September '03 new results have been 
published by the SNO Collaboration^, obtained after adding salt to their heavy wa- 
ter detector in order to increase the sensitivity to the neutral current channels. The 
previous results have been confirmed with increased accuracy. The allowed region for 
the LA solution has been further restricted with the elimination of the upper region 
in Am 2 2 . Of great importance have also been the first results from WMAP^ on the 
cosmic radiation background. The related determination of cosmological parameters, 
in combination with other measurements, leads to an upper limit on the cosmological 
neutrino density Vt v ^ 0.015. This is a very important result that indicates that 
neutrinos are not a major component of the dark matter in the Universe. For three 
degenerate neutrinos the WMAP limit implies an upper bound on the common mass 
given by m v < 0.23 eV Sl! . Given the priors that are assumed for this determination, 
i.e. a definite cosmological model, a 2-digit value for the bound is not to be taken 
too seriously EJ. Still the quoted value is about an order of magnitude smaller than 
the bound from tritium beta decay and of the same order of the upper bound on the 
Majorana mass that fixes the rate of neutrinoless double beta decay. 

In spite of this progress there are many alternative models of neutrino masses 
This variety is mostly due to the considerable experimental ambiguities that 
still exist. One first missing input is the absolute scale of neutrino masses: neu- 
trino oscillations only determine mass squared differences. For atmospheric neutrinos 
Am 2 tm ~ 2.6 10 -3 eV 2 while for solar neutrinos Am^, ~ 7 10~ 5 eV 2 . Another key 



missing quantity is the value of the third mixing angle S13 on which only a bound is 
known, s 13 < 0.22. Then it is essential to know whether the LSND signal which 
has not been confirmed by KARMEN and is currently being double-checked by 
MiniBoone , will be confirmed or will be excluded. If LSND is right we probably 
need at least four light neutrinos; if not we can do with only the three known ones. 

Here we will briefly summarize the main categories of neutrino mass models, dis- 
cuss their respective advantages and difficulties and give a number of examples. We 
illustrate how forthcoming experiments can discriminate among the various alterna- 
tives. We will devote a special attention to a comprehensive discussion in a GUT 
framework of neutrino masses together with all other fermion masses. This is for ex- 
ample possible in mo dels based on SU(5)x U(1)f or on SO(10) (we always consider 

susy gut's) Q3EE3. 



2. Basic Formulae and Data for Three-Neutrino Mixing 



We assume in the following that the LSND signal will not be confirmed so that 
there are only two distinct neutrino oscillation frequencies, the atmospheric and the 
solar frequencies. These two can be reproduced with the known three light neutrino 
species (for more than three neutrinos see, for example, ref. ^^). 

Neutrino oscillations are due to a misalignment between the flavour basis, v' = 
(i/ e , z/ M , u T ), where v e is the partner of the mass and flavour eigenstate e~ in a left- 
handed (LH) weak isospi n SU(2 ) doublet (similarly for and z/ r )) and the mass 
eigenstates v = (ui, u 2 , u 3 ) EHBJ : 



v 



where U is the unitary 3 by 3 mixing matrix. Given the definition of U and the 
transformation properties of the effective light neutrino mass matrix m u : 



iT 1 

U T m u U 



v T U T m v Uv 

Diag (mi, m 2 , m 3 ) = m diag 



(2) 



we obtain the general form of m v (i.e. of the light v mass matrix in the basis where 
the charged lepton mass is a diagonal matrix): 



m,. 



U*m diag U ] 



(3) 



The matrix U can be parameterized in terms of three mixing angles (9 12 , #23 and #13 
(0 < 9ij < 71"/ 2) and one phase <p (0 < <p < 2n) exactly as for the quark mixing 
matrix Vckm- The following definition of mixing angles can be adopted: 
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Table 1: Square mass differences and mixing angles. 
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5.4 


6.9 


9.5 


Aml tm (10- 3 eV 2 ) 


1.4 


2.6 


3.7 


sin 2 9i2 


0.23 


0.30 


0.39 


sin 2 9 2 3 


0.31 


0.52 


0.72 


sin 2 6>i3 





0.006 


0.054 



where = sin%, = cos%. In addition, if z/ are Majorana particles, we have the 
relative phases among the Majorana masses m 1; m 2 and ma. If we choose m 3 real 
and positive, these phases are carried by = Irai^e*^ 1,2 Thus, in general, 9 
parameters are added to the SM when non-vanishing neutrino masses are included: 
3 eigenvalues, 3 mixing angles and 3 CP violating phases. 

In our notation the two frequencies, Amj/4E (I = sun,atm), are parametrized 
in terms of the v mass eigenvalues by 

Am Ln = \ Am2 n\, Am ltm = \ Am l 3 \ ■ ( 5 ) 

where Am 2 2 = |m 2 | 2 — Im^ 2 > and Am^ = m\ — |m 2 | 2 . The numbering 1,2,3 
corresponds to our definition of the frequencies and in principle may not coincide with 
the ordering from the lightest to the heaviest state. From experiment, see table CpS, 
we know that S13 is small, according to CHOOZ, s 13 < 0.22 (3cr) ^D. Atmospheric 

neutrino oscillations mainly depend on (Am 2 tm , #23, #13), while solar oscillations are 
controlled by (Am 2 oi , 6*12, 6*13). Therefore, in the ideal limit of exactly vanishing S13, 
the solar and atmospheric oscillations decouple and depend on two separate sets 
of two-flavour parameters. For atmospheric neutrinos we have C23 ~ s 2 3 ~ 1/ \/2, 
corresponding to nearly maximal mixing. Oscillations of muon neutrinos into tau 
neutrinos are favoured over oscillations into sterile neutrinos (za,). The conversion 
probability and the zenith angular distribution of high-energy muon neutrinos are 



sensitive to matter effects, which distinguish v T from v s . Moreover, for conversion 
of into pure u s , neutral current events would become up/down asymmetric. In 
both cases data strongly disfavour the pure sterile case. Oscillations into u T are also 
indirectly supported by a SK data sample that can be interpreted in terms of enriched 
r-like charged-current events. The sterile component of the neutrino participating in 
atmospheric oscillations should amount to less than 0.25 at 90% C.L. Disappearance 
of laboratory-produced muon neutrinos has also been confirmed within expectations 
by the K2K experiment. 

The only surviving solution to the solar neutrino problem after KamLAND and 
SNO-salt results is LA' 22 ' 23 ' 2 ^ , with Am 2 sol « 7- 10" 5 eV 2 and sin 2 12 « 0.3. Before 
KamLAND the interpretation of solar neutrino data in terms of oscillations required 
the knowledge of the Boron neutrino flux, For instance, charged and neutral 
current data from SNO are sensitive, respectively, to fsiPee) an d to fB(J2 a Pea) (a = 
e,v,r), where (P e f) denotes the appropriately averaged conversion probability from 
v e to Uf. KamLAND^ provides a direct measurement of (P ee ). Beyond the impact 
on the oscillation parameters and a check that the solar standard model works well 
(/b = (1.00 ±0.06) x 5.05- 10 6 cm _2 sec _1 ), the comparison among these experiments 
shows that the conversion of Boron solar neutrinos into sterile neutrinos is compatible 
with zero. Therefore, the LSND indication for a third oscillation frequency associated 
to one or more sterile neutrinos is not supported by any other experiment, at the 
moment. Now, after KamLAND, also the possibility that such a frequency originates 
from a CPT violating neutrino spectrum has no independent support. Data from 
solar neutrino experiments and KamLAND, involving, respectively, electron neutrinos 
and electron antineutrinos, are compatible with a CPT invariant spectrum. 

If we take maximal S23 and keep only linear terms in u = Si3e lip , from experiment 
we find the following structure of the Ufi (/ = e,/x,r, i = 1, 2, 3) mixing matrix, apart 
from sign convention redefinitions: 

/ C12 _ S12 u \ 

U fi = \-(s 12 + c 12 u*)/V2 (C12 - s 12 u*)/y/2 1/V2 , (6) 
V (sia - Ci 2 «*)/V2 -(c 12 + s 12 u*)/V2 1/V2J 

where 6*12 is close to it/6 (for si 2 = l/v3 and u = we have the so-called tri- 
bimaximal mixing pattern ESJ, with the entries in the second column all equal to 
l/y/3 in absolute value). Given the observed frequencies and our notation in eq. (jSJ, 
there are three possible patterns of mass eigenvalues: 



Degenerate 
Inverted hierarchy 
Normal hierarchy 



|mi| ~ \m 2 \ ~ |m 3 | ^> \rrii — nij\ 
\mi\ ~ \m 2 \ 3> |m3| 

\m 3 \ > |m 2 ,i| (7) 



Models based on all these patterns have been proposed and studied and all are in 
fact viable at present. In the following we will first discuss neutrino masses in general 



and, in particular, Majorana neutrinos. Then we recall the existing constraints on 
the absolute scale of neutrino masses. We then discuss the importance of neutrinoless 
double beta decay that, if observed, would confirm the Majorana nature of neutrinos. 
Also the knowledge of the rate of this process could discriminate among the possi- 
ble patterns of neutrino masses in (JJJ). The possible importance of heavy Majorana 
neutrinos for the explanation of baryogenesis through leptogenesis in the early Uni- 
verse will be briefly discussed. We finally review the phenomenology of neutrino mass 
models based on the three spectral patterns in (J7J) and the respective advantages and 
problems. 



3. Neutrino Masses and Lepton Number Violation 

Neutrino oscillations imply neutrino masses which in turn demand either the ex- 
istence of right-handed (RH) neutrinos (Dirac masses) or lepton number L violation 
(Majorana masses) or both. Given that neutrino masses are certainly extremely 
small, it is really difficult from the theory point of view to avoid the conclusion that 
L conservation must be violated. In fact, in terms of lepton number violation the 
smallness of neutrino masses can be explained as inversely proportional to the very 
large scale where L is violated, of order Mqut or even Mp\. 

Once we accept L non- conservation we gain an elegant explanation for the small- 
ness of neutrino masses. If L is not conserved, even in the absence of heavy RH 
neutrinos, Majorana masses for neutrinos can be generated by dimension five opera- 
tors f th e form 

0s = mME!k +hc (8) 

with H being the ordinary Higgs doublet, li the SU(2) lepton doublets, A a matrix in 
flavour space, A a large scale of mass, of order Mqut or Mpi and a charge conjugation 
matrix C between the lepton fields is understood. Neutrino masses generated by O5 
are of the order m v « v 2 /A for A^ ~ 0(1), where v ~ 0(100 GeV) is the vacuum 
expectation value of the ordinary Higgs. 

We consider that the existence of RH neutrinos v c is quite plausible because all 
GUT groups larger than SU(5) require them. In particular the fact that u c completes 
the representation 16 of SO(10): 16=5+10+1, so that all fermions of each family are 
contained in a single representation of the unifying group, is too impressive not to 
be significant. At least as a classification group SO(10) must be of some relevance. 
Thus in the following we assume that there are both v c and L non- conservation. With 
these assumptions the see-saw mechanism^ is possible. Also to fix notations we recall 
that in its simplest form it arises as follows. Consider the SU(3) x SU(2) x U(l) 
invariant Lagrangian giving rise to Dirac and v c Majorana masses (for the time being 
we consider the v (versus v c ) Majorana mass terms as comparatively negligible): 

£ = -y cT y u (Hl) + l -v cT Mv c + h.c. (9) 



The Dirac mass matrix uid = Vuv/y/2, originating from electroweak symmetry break- 
ing, is, in general, non-hermitian and non- symmetric, while the Majorana mass matrix 
M is symmetric, M = M T . We expect the eigenvalues of M to be of order Mqut 
or more because v c Majorana masses are SU(3)x SU(2)x U(l) invariant, hence un- 
protected and naturally of the order of the cutoff of the low-energy theory. Since all 
v c are very heavy we can integrate them away. For this purpose we write down the 
equations of motion for v c in the static limit, i.e. neglecting their kinetic terms: 

dC 

- — = y u (Hi)-Mv c = . (10) 



From this, by solving for u c , we obtain: 



M~ x y v [EV) . (11) 



We now replace in the lagrangian, eq. Q, this expression for v c and we get the 
operator O5 of eq. (jHJ) with 

y = -y, T M-V , (12) 

and the resulting neutrino mass matrix reads: 

m v = m T D M- l m D . (13) 

This is the well known see-saw mechanism result ^ : the light neutrino masses are 
quadratic in the Dirac masses and inversely proportional to the large Majorana mass. 
If some v c are massless or light they would not be integrated away but simply added 
to the light neutrinos. Notice that the above results hold true for any number n of 
heavy neutral fermions R coupled to the 3 known neutrinos. In this more general case 
M is an n by n symmetric matrix and the coupling between heavy and light fields is 



described by the rectangular n by 3 matrix m D . Note that for m v « J Am^ tm ~ 0.05 
eV and m u « m 2 D /M with mo ~ v « 200 GeV we find M w 10 15 GeV which indeed 
is an impressive indication for Mqut- 

If additional non-renormalizable contributions to O5, eq. (jHJ), are comparatively 
non-negligible, they should simply be added. For instance in SO(10) or in left-right 
extensions of the SM, an SU(2) l triplet can couple to lepton doublets and may induce 
a sizeable contribution to neutrino masses. At the level of the low-energy effective 
theory, such contribution is still described by the operator O5 of eq. (JHJ), obtained 
by integrating out the heavy SU(2)l triplet (which also acquires a VEV due to its 
coupling to the Higgs doublets). This contribution is called type II to be distinguished 
from that obtained by the exchange of RH neutrinos (type I). After elimination of 
the heavy fields, at the level of the effective low-energy theory, the two types of terms 
are equivalent. In particular they have identical transformation properties under a 



chiral change of basis in flavour space. The difference is, however, that in the see-saw 
mechanism, the Dirac matrix mo is presumably related to ordinary fermion masses 
because they are both generated by the Higgs mechanism and both must obey GUT- 
induced constraints. Thus if we assume the see-saw mechanism in its simplest type I 
version, more constraints are implied. 

4. Importance of Neutrinoless Double Beta Decay 

Oscillation experiments do not provide information about the absolute neutrino 
spectrum and cannot distinguish between pure Dirac and Majorana neutrinos. From 
the endpoint of tritium beta decay spectrum we have an absolute upper limit of 
2.2 eV (at 95% C.L.) on the mass of electron antineutrino , which, combined with 
the observed oscillation frequencies under the assumption of three CPT-invariant light 
neutrinos, represents also an upper bound on the masses of the other active neutrinos. 
Complementary information on the sum of neutrino masses is also provided by the 
galaxy power spectrum combined with measurements of the cosmic microwave back- 
ground anisotropics. According to the recent analysis of the WMAP collaboration^, 
J2i \ m i\ < 0.69 eV (at 95% C.L.). More conservative analyses^ give J2i \ m i\ < 1-01 
eV, still much more restrictive than the laboratory bound. 

The discovery of 0vj3(3 decay would be very important because it would establish 
lepton number violation and the Majorana nature of u's, and provide direct informa- 
tion on the absolute scale of neutrino masses. As already mentioned the present limit 
from QvfiP is \m ee \ < 0.2 eV or to be more conservative \m ee \ < (0.3^-0.5) eV 1 29 1 30 } 
Note, however, that the WMAP limit implies for 3 degenerate i/'s \m\ < 0.23 eV and, 
taken at face value, this limit would pose a direct constraint on m ee . 

It is interesting to see what is the level at which a si gnal c an be expected or at 
least not excluded in the different classes of models in Q 131 1 32] _ q Uan tity which 
is bound by experiments is the 11 entry of the v mass matrix, which in general, from 
eqs. (J2J) and (HJ), is given by : 

\m ee \ = |(1 - s\z) (raic^ + m 2 s 2 12 ) + m 3 e 2l<l> s 2 13 \ , (14) 

For 3-neutrino models with degenerate, inverse hierarchy or normal hierarchy mass 
patterns, starting from this general formula it is simple to derive the following bounds. 

a) Degenerate case. If |m| is the common mass and we take S13 = 0, which is a safe 
approximation in this case, because \m 3 \ cannot compensate for the smallness 
of S13, we have m ee ~ \m\(c\ 2 ±s\ 2 ). Here the phase ambiguity has been reduced 
to a sign ambiguity which is sufficient for deriving bounds. So, depending on 
the sign we have m ee = \m\ or m ee = \m\cos26i 2 . We conclude that in this 
case m ee could be as large as the present experimental limit but should be at 
least of order 0(yAm^ ra ) ~ O(10~ 2 eV) unless the solar angle is practically 



maximal, in which case the minus sign option can be arbitrarily small. But 
the experimental 2-cr range of the solar angle does not favour a cancellation by 
more than a factor of 3. 

b) Inverse hierarchy case. In this case the same approximate formula m ee = 
|m|(c^ 2 ±s^ 2 ) holds because m 3 is small and S13 can be neglected. The difference 
is that here we know that \m\ ~ J lS.vn? atm so that \m ee \ < \J Am„ ta ~ 0.05 
eV. At the same time, since a full cancellation between the two contributions 
cannot take place, we expect \m ee \ > 0.01 eV. 

c) Normal hierarchy case. Here we cannot in general neglect the m.3 term. However 
in this case \m ee \ ~ | Am 2 sun s\ 2 ± \J Aml tm sf 3 | and we have the bound 
\m ee \ < a few 10~ 3 eV. 

Recently evidence for 0^/3/3 was claimed in ref. ^3 a ^ the 4.2-cx level corresponding 
to \m ee \ ~ (0.2-=-0.6) eV ((O.l-hO.9) eV in a more conservative estimate of the involved 
nuclear matrix elements). If confirmed this would rule out cases b) and c) and point to 
case a) or to models with more than 3 neutrinos. We recall that further contributions 
to 0^/3/3 transition amplitudes might occur in models with additional L violating 
interactions, such as R-parity breaking supersymmetry. 

5. Baryogenesis via Leptogenesis from Heavy v c Decay 

In the Universe we observe an apparent excess of baryons over antibaryons. It is 
appealing that one can explain the observed baryon asymmetry by dynamical evo- 
lution (baryogenesis) starting from an initial state of the Universe with zero baryon 
number. For baryogenesis one needs the three famous Sakharov conditions: B viola- 
tion, CP violation and no thermal equilibrium. In the history of the Universe these 
necessary requirements can have occurred at different epochs. Note however that 
the asymmetry generated by one epoch could be erased at following epochs if not 
protected by some dynamical reason. In principle these conditions could be verified 
in the SM at the electroweak phase transition. B is violated by instantons when kT 
is of the order of the weak scale (but B-L is conserved), CP is violated by the CKM 
phase and sufficiently marked out-of- equilibrium conditions could be realized during 
the electroweak phase transition. So the conditions for baryogenesis at the weak scale 
in the SM superficially appear to be present. However, a more quantitative analysis 

shows that baryogenesis is not possible in the SM because there is not enough 
CP violation and the phase transition is not sufficiently strong first order, unless 
van < 80 GeV, which is by now completely excluded by LEP. In SUSY extensions 
of the SM, in particular in the MSSM, there are additional sources of CP violation 
and the bound on rrtu is modified by a sufficient amount by the presence of scalars 
with large couplings to the Higgs sector, typically the s-top. What is required is that 



rrih ~ 80 — 110 GeV, a s-top not heavier than the top quark and, preferentially, a 
small tan/3. But also this possibility has by now become at best marginal with the 
results from LEP2. 

If baryogenesis at the weak scale is excluded by the data it can occur at or just 
below the GUT scale, after inflation. But only that part with |B — L| > would 
survive and not be erased at the weak scale by instanton effects. Thus baryogenesis 
at kT ~ 10 10 — 10 15 GeV needs B-L violation at some stage like for m v if neutrinos 
are Majorana particles. The two effects could be related if baryogenesis arises from 
leptogenesis then converted into baryogenesis by instantons Recent results on 
neutrino masses are compatible with this elegant possibility^^. Thus the case of 
baryogenesis through leptogenesis has been boosted by the recent results on neutrinos 

ma 

In leptogenesis the baryon asymmetry is produced by the out of equilibrium, CP 
and L-violating decays of heavy right-handed neutrinos v c . In the simplest cases the 
spectrum of RH neutrinos is assumed to be hierarchical, M\ -C M.2$ and the mecha- 
nism is dominated by the lightest state, v\. Thus the expected baryon asymmetry is 
proportional to the product e^b between the CP decay asymmetry 

ei -rK-/) + rK-/) ' (15) 

and the efficiency factor 8 < 1, the fraction of the produced asymmetry that survives 
after v\ decay. To keep S close to 1 and avoid washing out the developed lepton asym- 
metry, L violating interactions should be sufficiently weak to stay out of equilibrium 
when vl decays. This in turns happens when the lifetime of v\ exceeds the age of the 
Universe, which is expressed by the condition 

m x = {mD t D)lX < ID" 3 eV , (16) 
Mi y J 

with min{|ml|, \m2\, |m3|} < mj. This condition does not provide yet an abso- 
lute bound on neutrino masses since, in realistic simulations, the observed baryon 
asymmetry is achieved for S less than one. For hierarchical RH neutrino masses the 
asymmetry ei is given by: 

e l *- 3 £ T J 1 -j:irn[(y„yi)l j \^ , (17) 

(a = 1 in the supersymmetric case and 1/2 in the non-supersymmetric one) which 
gives rise to the Davidson-Ibarra (DI) bound 

lei I < — . * (max{|ml|, \m2\, \m3\} — min{|ml|, \m2\, \m3\}) , (18) 
Stt (-n ) 



where (H°) denotes the VEV of the Higgs doublet giving mass to the up quarks. 
Several interesting constraints on the neutrino spectrum emerge from the above dis- 
cussion. First of all, in the assumed limit 3> Ml, the CP asymmetry vanishes 
for a completely degenerate light neutrino spectrum (of course this limit is rather 
unnatural in the see-saw context). Moreover, depending on the details of the as- 
sumed cosmological scenario which fixes the initial abundance of RH neutrinos, the 
DI bound translates into a lower bound o n the l ightest RH neutrino mass Mi, ranging 
from approximately 10 7 GeV to 10 9 GeV^^^. Finally, an absolute upper bound on 
light neutrino masses can be derived. Indeed, if we enhance the absolute scale of light 
neutrino masses, on the one hand we depart from the out-of-equilibrium condition of 
eq. (fTHjl and on the other hand the DI bound, proportional to Am^ tm /m 3 , becomes 
more and more stringent. Quantitative studies'"^ show that 

|m;| < (0.12 -=-0.15) eV . (19) 

It should be stressed that these bounds hold under the assumption of strict hierarchy 
in the RH neutrino sector. Indeed, by relaxing the assumption Ms^ 3> Mi, there are 
important corrections to the expression of e\ in eq. (fTTj) . which may receive a resonant 
enhancement and which do not vanish any longer for degenerate light neutrinos. For 
a moderate degeneracy among RH neutrinos, the DI bound is violated and Mi can 
be considerably lower than 10 7 — 10 9 GeV^*. At the same time, the enhanced CP 
violating asymmetry allows a successful leptogenesis even for light neutrino masses 
around the eV scale In any case it is impressive that the resulting range of 

neutrino masses is fully consistent with the results on neutrino oscillations. 



6. Degenerate Neutrinos 

For degenerate neutrinos the average m 2 is much larger than the splittings. At first 
sight the degenerate case is the most appealing: the observation of nearly maximal 
atmospheric neutrino mixing and the more recent result that also the solar mixing 
is large suggests that all v masses are nearly degenerate. We shall see that this 
possibility has become less attractive with the recent new experimental information. 

It is clear that in the degenerate case the most likely origin of v masses is from 
some dimension 5 operators (Hl)J Xij(Hl) j j A not related to the see-saw mechanism 
m„ = m^)M~ l m£)- In fact we expect the v Dirac mass mo not to be degenerate like for 
all other fermions and a conspiracy to reinstate a nearly perfect degeneracy between 
m D and M, which arise from completely different physics, looks very unplausible (see, 
however, ED). Thus in degenerate models, in general, there is no direct relation with 
Dirac masses of quarks and leptons and the possibility of a simultaneous description 
of all fermion masses within a grand unified theory is more remote l 4 ^ 4 *^ . 

The degeneracy of neutrinos should be guaranteed by some slightly broken sym- 
metry. Models based on discrete or continuous symmetries have been proposed. For 



example in the models of ref. ' ' ' the symmetry is SO (3): in the unbroken limit 
neutrinos are degenerate and charged leptons are massless. When the symmetry is 
broken the charged lepton masses are much larger than neutrino splittings because 
the former are first order while the latter are second order in the electroweak symme- 
try breaking. In this kind of models the mixing angles are completely undetermined 
in the symmetric phase and they originate only in the spontaneously broken phase 
from a misalignment between the symmetry breaking terms for neutrinos and charged 
leptons. 

In principle, when considering models with degenerate masses we must keep in 
mind that radiative corrections can modify mass splittings and mixing angles in the 
running from the high scale where neutrino masses are determined at the fundamental 
level (i.e. the heavy Majorana mass M or Mqut) down to the electroweak scale 
These running effects can be evaluated by renormalisation group techniques. The 
effects depend on the parameters A ab = (m a + m h ) / (m a — m h ) and are, with good 
accuracy, proportional to e = ?/ 2 /(167r 2 ) log(M/m^), where y T is the r-lepton Yukawa 
coupling. The value of e is around 10~ 5 in the SM and larger by a factor 1 + tan 2 j3 
in the MSSM. The corrections are negligible for \A a t,e\ ^ 0(1). When some of these 
quantities are large a rapid transition takes place towards a fixed point configuration 
of mixing angles that does not correspond to the observed pattern. In practice, given 
the present upper bounds on the degenerate neutrino common mass mo and the LA 
value of Am 2 oi , the corrections are always negligible in the SM and can only become 
sizable in the MSSM if tan 2 f3 is large, mo is close to its absolute upper bound and m\ 
and mi (those entering in the smallest difference Am 2 oZ ) are nearly equal in absolute 
value and sign. Note that these remarks do not include the effects from thresholds 
that could affect running in a significant way, depending on the details of the heavy 
and/or light particle spectrum. 

The upper limit on the common value \m\ becomes particularly stringent if one 
adopts the cosmological WMAP bound |m| < 0.23 eV ^' (or the more conservative 
one \m\ < 0.34 eV^). The more direct laboratory limit from tritium beta decay 
is \m\ < 2.2 eV In past years degenerate models with v masses as large as 

\m\ ~ (1 -T- 2) eV were considered with the perspective of a large fraction of hot dark 
matter in the universe. In this case, however, the existing limit on the abs ence of 
OvfiP (|m ee | < 0.2 eV or to be more conservative \m ee \ < (0.3-^0.5) eV) implies E7 | 30 ) 
approximate double maximal mixing (bimixing) for solar and atmospheric neutrinos. 
As discussed in sect. 4 , for |m| 3> m ee , one needs mi ~ — m 2 and, to a good accuracy, 
c i2 ~ s i2' m order to satisfy the bound on m ee . This is exemplified by the following 
texture 

/ -l/y/2 l/y/2 \ 

m u = m\ -l/y/2 (l + r/)/2 (1 + ^/2 , (20) 
V l/y/2 (l + tj)/2 {1 + V)/2J 

where r|<Cl, corresponding to an exact bimaximal mixing, S13 = and the eigenval- 



ues are mi = m, 7712 = —m and 772,3 — (1+77)771. This texture has been proposed in the 
context of a spontaneously broken SO (3) flavor symmetry and it has been stu died to 
analyze the stability of the degenerate spectrum against radiative corrections H^ l ^bj _ 
A more realistic mass matrix can be obtained by adding small perturbations to m v 
in eq. Ij21ty : 



771,, = 771 



/ 5 -1/V2 (l-e)/V2 

-1/V2 (1 + ?7)/2 (l + V -e)/2 1 (21) 
V (l-e)/V2 (l + 7/-e)/2 (l + 7j-2e)/2 



where e parametrizes the leading flavor- dependent radiative corrections (mainly in- 
duced by the r Yukawa coupling) and 5 controls m ee . Consider first the case 5 <C e. 
To first approximation # 12 remains maximal. We get Am? sun m m 2 e 2 /rj and 

~ {a^) • m " m { S5 J • (22) 

If we instead assume 5 > e, we find &m 2 sun ss 2m 2 5, 2 3 ~ tt/4, sin 2 20 12 ~ 1 — <5 2 /4. 
Also in this case the solar mixing angle remains unacceptably close to 7r/4, unless 
further contributions to the mixing matrix are induced from the diagonalization of 
the charged lepton sector. We get: 

Atti 2 

13 «O , T71 ee ^— , (23) 

2771 

too small for detection if the average neutrino mass m is around the eV scale. We see 
that with increasing \m\ more and more fine tuning is needed to reproduce the LA 
solution values of A77i 2 un and 612 ■ In conclusion, even without invoking the WMAP 
limit, large v masses, |m| ~ (1^-2) eV, are disfavoured by the limit on 0v/3[3 decay 
and by the emerging of the LA solution with the solar angle definitely not maximal. 
From these considerations it is possible to derive the bound \m\ < 0.9 h eV (90% 
C.L.) where h ~ 1 parametrizes the uncertainties in the nuclear matrix elements. 
Also, we have seen in sect. 5 that the attractive mechanism of baryogenesis through 
leptogenesis appears to disfavour \m\ ^ 0.1 eV, at least in the simplest realizations. 
All together, after WMAP and KamLAND, among degenerate models those with 
|m| ^ (0.23 -T- 1) eV are favoured by converging evidence from different points of 
view. 

For I m I not larger than the 0^/3/3 bound one does not need a cancellation in m ee 
and mi and 7712 can be approximately equal in magnitude and phase. For example, 
in the limit S13 = 0, the matrix 



m v = m -1 , (24) 




corresponds to maximal #23 (pseudo Dirac 23 sub- matrix) with Diag[m v ] — m(l, 1, — 1) 
The angle 9 12 is unstable and a small perturbation can give any value to it. Note, 
however, that in this case the non vanishing matrix elements must be of equal abso- 
lute value and not just of order 1. So either this is guaranteed by a symmetry (as, 
for example, in ref. or the model is unnaturally fine-tuned. 

As a different example (also with no cancellation between m\ and mg), a model, 
which is simple to describe but difficult to derive in a natural way, is one I ^ I ^Q where 
up quarks, down quarks and charged leptons have "democratic" mass matrices, with 
all entries equal (in first approximation): 



rrii 




+ 5m f , (25) 



where rhj (f = u,d,e) are three overall mass parameters and 5mj denote small 
perturbations. If we neglect 5m f, the eigenvalues of mj are given by (0, 0,3 m/). The 
mass matrix mj is diagonalized by a unitary matrix Uf which is in part determined 
by the small term 5m f. If 5m u « 5md, the CKM matrix, given by Vckm — U^Ud, is 
nearly diagonal, due to a compensation between the large mixings contained in U u and 
Ud- When the small terms 8m f are diagonal and of the form 5m f = Diag(— e/, e/, 8f) 
with 5f ^> 6f, the matrices Uf are approximately given by (note the analogy with the 
quark model eigenvalues 7r°, 77 and 77'): 



(l/y/2 -l/y/2 

1/V6 1/y/E -2/y/E I . (26) 
V1/V3 1/V3 i/Vs 



Note that, due to the degeneracy in the 1,2 sector in the unperturbed limit, any 
superposition of the first two rows would also be an eigenvector of zero mass. Thus 
the choice (C//)i3 = is unjustified in the unperturbed limit and is only determined 
by a particular form of the perturbation. At the same time, the lightest quarks and 
charged leptons acquire a non-vanishing mass. The leading part of the mass matrix 
in eq. (|23j) is invariant under a discrete S 3 l x permutation symmetry. The same 
requirement leads to the general neutrino mass matrix: 



m,. = m 





+ 5m v , (27) 



where 5m v is a small symmetry breaking term and the two independent invariants 
are allowed by the Majorana nature of the light neutrinos. If r vanishes the neutrinos 
are almost degenerate. In the presence of 5m v the permutation symmetry is broken 
and the degeneracy is removed. If, for example, we choose 5m v = Diag(0, e, rj) with 
e < 77 1 and r e, the solar and the atmospheric oscillation frequencies are 



determined by e and 77, respectively. The problems with this class of models are 
that the solar angle should be close to maximal (typically more than the atmospheric 
angle) and that the neutrino spectrum and mixing angles are not determined by the 
symmetric limit (this applies, for example, to S13 ~ 0) but only by a specific choice of 
the parameter r and of the perturbations that cannot be easily justified on theoretical 
grounds. Notice also that the simplest choice of parameters leads to sin 2 26* 2 3 very 
close to 8/9, value which is now excluded at the 2 a level. 

In this model the mixing angles are almost entirely due to the charged lepton 
sector. Recently the question of whether observed neutrino mixings can dominantly 
arise from the charged lepton sector in a natural way was studied in general in ref. 

Of course, one can always choose an ad hoc basis where this is true: the point 
is to decide whether this formal choice can be naturally justified in the physical 
basis where the symmetries of the lagrangian are specified. The conclusion is that 
in presence of two large mixing angles 812 and 6*23 with the third angle #13 being 
small, the construction of a natural model with dominance of U e is made much more 
difficult than in the case of only t he atmosph eric angle 823 large. Exemples of natural 
models of this sort can be given | 45|50 | 51 |52j anc j stated difficulty is reflected in 
the relatively complicated symmetry structure required. 

In conclusion, the parameter space for degenerate models has recently become 
smaller because of the indications from WMAP (and also, too some extent, from lep- 
togenesis) that tend to lower the maximum common mass allowed for light neutrinos. 
It is also rather difficult to reproduce the observed pattern of frequencies and mix- 
ing angles, in particular two large and one small mixing angle, with the solar angle 
large but not maximal. Degenerate models that fit can only arise from a very special 
dynamics or a non abelian flavour symmetry with suitable breakings. 



6.1. Anarchy 

Anarchical models can be considered as particular cases of degenerate models 
with m 2 ~ Am 2 tm . In this class of models mass degeneracy is replaced by the principle 
that all mass matrices are structure-less in the neutrino sector, (including the LH 
charged fermions and possibly the RH neutrinos). For the LA solution the ratio of 
the solar and atmospheric frequencies is not so small: r = (AT7i 2 un )^^/A?7i 2 tm ~ 1/40 
and two out of three mixing angles are large. The key observation is that the see-saw 
mechanism tends to enhance the ratio of eigenvalues: it is quadratic in so that a 
hierarchy factor / in mc becomes f 2 in m u and the presence of the Majorana matrix 
M results in a further widening of the distribution. Another squaring takes place in 
going from the masses to the oscillation frequencies which are quadratic. As a result, 
a random generation of the and M matrix elements leads to a distribution of 
r that peaks around 0.1. At the same time the distribution of sin 2 8ij is rather flat 
for all three mixing angles. Clearly the smallness of #13 is problematic for anarchy. 



This can be turned into the prediction that in anarchical models #13 must be near the 
present bound (after all the value 0.2 for sin #13 is not that smaller than the maximal 
value 0.707). In conclusion, if 6*13 is near the present bound, one can argue that the 
neutrino masses and mixings, interpreted by the see-saw mechanism, can just arise 
from structure-less underlying Dirac and Major ana matrices. 



7. Inverted Hierarchy 



The inverted hierarchy configuration \m,\\ ~ | rri2 1 S> |tt^3 | consists of two levels 
mi and m-i with small splitting Am^ = /Am 2 sun and a common mass given by 
\m\ 2 \ ~ \/S.m 2 atm \ ~ 2.6 • 10~ 3 eV 2 . One particularly interesting example of this 
sort^*, which leads to double maximal mixing, is obtained with the phase choice 
m\ = —vri2 so that, approximately: 



m d i a g = Diag(\/2m, -V2m, 0) . (28) 

The effective light neutrino mass matrix 

m u = U*m diag U^ , (29) 

which corresponds to the mixing matrix of double maximal mixing C12 = S12 = 1/ 
and S13 = u = in eq. (JBJ). 

[1/V2 l/y/2 \ 

U fi = -1/2 1/2 l/y/2 , (30) 
V 1/2 -1/2 l/y/2 J 




m v = m -1 . (31) 



is given by: 



This texture for m v can be reproduced by imposing a U(l) flavour symmetry with 
charge L e — — L T starting either from (Hl)J Xij(Hl) j / A or from RH neutrinos via 
the see-saw mechanism. However the absolute values of the 12 and 13 terms would be 
in general different in this case. As a consequence, while the vanishing of S13 and the 
maximal value of 9 12 are still valid, the atmospheric angle deviates from the maximal 
value with tan #23 = x where x is the absolute value of the ratio Tn ul3 /m ul o. We also 
note that the 1 — 2 degeneracy remains stable under radiative corrections HHESl_ 
The leading texture in (J3T]) can be perturbed by adding small terms: 



m u = m — 1 77 77 , (32) 




where 5 and rj are small (<^C 1), real parameters defined up to coefficients of order 
1 that can differ in the various matrix elements. One could also make the absolute 
values of the 12, 13 terms different by terms of order S\s5. The perturbations leave 
Am 2 tm and #23 unchanged, in first approximation. We obtain tan 2 9 12 ~ 1 + S + rj and 
Am 2 sun / Am 2 atm »s i] + 5, where coefficients of order one have been neglected. Moreover 
9i 3 ~ 77. If r] 3> 5, we have 



Ain't. 



'13 



In the other case, 77 <C 5 we obtain: 



013 < 



Am 2 



atm 



Am 2 . 



Am 2 



atm , 



Am2 sun 



y Am 2 tm ^ 



(33) 



(34) 



Indeed, 
close to 



There is a well-known difficulty of this scenario to fit the LA solution ^^ 1 ^ 
barring cancellation between the perturbations, in order to obtain a Am 2 nr 
the best fit LA value, r? and 5 should be smaller than about 0.1 and this keeps the 
value of sin 2 2^12 very close to 1, in disagreement with global fits of solar data^T 
Notice that the required deviation of the solar mixing angle from the maximal value is 
of the order of the Cabibbo angle 9c and indeed the empirical relation 9i 2 + 9 c = 7r/4 
holds within the experimental errors However, even starting from exact bimixing 
the pattern of parameters needed to bring the solar angle down from the maximal 
value can be obtained in a natural way as an effect of the charge d lepton matrix 
diagonalization . This possibility, studied in detail in refs. EDEU^ j s no ^ excluded 
but is strongly constrained by the observed smallness of S13, as in general the amount 
of deviation from maximal solar angle is typically of order S13 (while the deviation 
from maximal atmospheric mixing are of second order). 

With the phase choice mi = i. e. for Diag[m 
S13 = 0, one obtains the matrix 

/ 1 

m v = m 1/2 -1/2 
V -1/2 1/2 

which corresponds to large #23 with the solar angle unstable to small perturbations. 
In this case fine tuning or a symmetry is necessary to fix the ratios of the matrix 
elements as indicated. 

In conclusion, also for inverse hierarchy some special dynamics or symmetry is 
needed to reproduce in detail the observed features of the data. 



m(l,l,0), in the limit 



(35) 



8. Normal Hierarchy 



We now discuss the class of models which we consider the simplest approach to 
neutrino masses and mixings. In particular, in this context one can formulate the most 



constrained framework which allows a comprehensive combined study of all fermion 
masses in GUT's. We start by assuming three widely split v's and the existence of 
a RH neutrino for each generation, as required to complete a 16 dimensional repre- 
sentation of SO (10) for each generation. We then assume dominance of the see-saw 
mechanism m v = rn^M^m^. We know that the third-generation eigenvalue of the 
Dirac mass matrices of up and down quarks and of charged leptons is systematically 
the largest one. It is natural (although not necessary) to imagine that this property 
could also be true for the Dirac mass of z/s: m d ^ a9 ~ Diag(0, 0, moz)- After see-saw 
we expect m u to be even more hierarchical being quadratic in (barring fine- 
tuned compensations between mp and M). Note however that, for the LA solution: 
r ~ 1/40, so that the required amount of hierarchy, r = /\m 2 atm / /S.m 2 sun = ml/m 2 , is 
quite moderate. 

A possible difficulty for hierarchical models is that one is used to expect that 
large splittings correspond to small mixings because normally only close-by states are 
strongly mixed. In a 2 by 2 matrix context the requirement of large splitting and 
large mixings leads to a condition of vanishing determinant and large off-diagonal 
elements. For example the matrix 

(* l) < 36 > 

has eigenvalues and 1 + x 2 and for x of 0(1) the mixing is large. Thus in the 
limit of neglecting small mass terms of order m^2 the demands of large atmospheric 
neutrino mixing and dominance of m 3 translate into the condition that the 2 by 2 
sub determinant 23 of the 3 by 3 mixing matrix approximately vanishes. The problem 
is to show that this vanishing can be arranged in a natural way without fine tuning. 
Once near maximal atmospheric neutrino mixing is reproduced the solar neutrino 
mixing can be arranged to be either small or large without difficulty by implementing 
suitable relations among the small mass terms. 

It is not difficult to imagine mechanisms that na turally lead to the approximate 
vanishing of the 23 sub-determinant. For example inEOEIJ \i j s assumed that one v c 
is particularly light and coupled to fi and r. In a 2 by 2 simplified context if we have 

M K (J J) , M-^A I) , mj ,= £ J) , ,37) 

then for a generic nip we find 

m v = vdF D M- x m D « ~ {\ ^) • (38) 

A different possibility that we find attractive is that, in the limit of neglecting terms 
of order m 12 and, in the basis where charged leptons are diagonal, the Dirac matrix 



ran, defined by u c mr,v, takes the approximate form, called "lopsided" 




/0 0\ 
mj)« 
\0 x l) 



(39) 



This matrix has the property that for a generic Majorana matrix M one finds: 



The only condition on M _1 is that the 33 entry is non zero. However, when the 
approximately vanishing matrix elements are replaced by small terms, one must also 
assume that no new 0(1) terms are generated in m u by a compensation between small 
terms in and large terms in M _1 . 

It is important for the following discussion to observe that mo given by eq. (}3*9"j) 
under a change of basis transforms a.s mo —>■ V^tudU where V and U rotate the right 
and left fields respectively. It is easy to check that in order to make mo diagonal we 
need large left mixings (i.e. large off diagonal terms in the matrix that rotates LH 
fields). Thus the question is how to reconcile large LH mixings in the leptonic sector 
with the observed near diagonal form of Vckm, the quark mixing matrix. Strictly 
speaking, since Vckm = U^Ud, the individual matrices U u and Ud need not be near 
diagonal, but Vckm does, while the analogue for leptons apparently cannot be near 
diagonal. However for quarks nothing forbids that, in the basis where m u is diagonal, 
the d quark matrix has large non diagonal terms that can be rotated away by a pure 
RH rotation. We suggest that this is so and that in some way RH mixings for quarks 
correspond to LH mixings for leptons. 

In the context of (SUSY) SU(5) there is a very attractive hint of how this sort of 
mechanism can be realized EMD. In t h e 5 of SU(5) the d c singlet appears together 
with the lepton doublet (u,e). The (u,d) doublet and e c belong to the 10 and v c to 
the 1 and similarly for the other families. As a consequence, in the simplest model 
with mass terms arising from only Higgs pentaplets, the Dirac matrix of down quarks 
is the transpose of the charged lepton matrix: rrid = {rrii) T . Thus, indeed, a large 
mixing for RH down quarks corresponds to a large LH mixing for charged leptons. 
At leading order we may have the lopsided texture: 




(40) 




/0 
1 
\0 1 




(41) 



In the same simplest approximation with 5 or 5 Higgs, the up quark mass matrix 
is symmetric, so that left and right mixing matrices are equal in this case. Then 
small mixings for up quarks and small LH mixings for down quarks are sufficient to 



guarantee small Vqkm mixing angles even for large d quark RH mixings. It is well 
known that a model where the down and the charged lepton matrices are exactly the 
transpose of one another cannot be exactly true because of the e/d and fi/s mass 
ratios. It is also known that one remedy to this problem is to add some Higgs com- 
ponent in the 45 representation of SU(5) But the symmetry under transposition 
can still be a good guideline if we are only interested in the order of magnitude of the 
matrix entries and not in their exact values. In models with U(l)p flavour symmetry 
equal second and third generation charges for the 5 could induce a lopsided form for 
both the charged leptons and the Dirac neutrino mass matrices. Otherwise, in the 
very crude model where the Higgs pentaplets come from a pure 10 representation of 
SO (10) one has nip = m u , i. e. the Dirac neutrino mass matrix is the same as 
the up quark mass matrix . For mu the dominance of the third family eigenvalue as 
well as a near diagonal form could be an order of magnitude remnant of this broken 
symmetry. 

To get a realistic mass matrix, we allow for deviations from the symmetric limit 
of ()40j) with x ~ o(l). For instance, we can consider those models where the neutrino 
mass matrix elements are dominated, via the see-saw mechanism, by the exchange of 
two right-handed neutrinos'^. Since the exchange of a single RH neutrino gives a 
successful zeroth order texture, we are encouraged to continue along this line. Thus, 
we add a sub-dominant contribution of a second RH neutrino, assuming that the third 
one gives a negligible contribution to the neutrino mass matrix, because it has much 
smaller Yukawa couplings or is much heavier than the first two. The Lagrangian that 
describes this plausible subset of see-saw models, written in the mass eigenstate basis 
of RH neutrinos and charged leptons, is 

C = y i u c Hl l + yy / Hh + ^-u c2 + ^-u c ' 2 , (42) 



/ \ ViVj . ViVj /ao\ 

\m v )ij oc — - H , 43 

v uj% 3 M M' y J 



leading to 

where i,j = {e, /i, r}. In particular, if y e <C ps y T and y'^ ~ y' T , we obtain: 

/ 5 e e \ 

m u = m e 1+7] 1 + r] \ , (44) 




where coefficients of order one multiplying the small quantities 5, e and r\ have been 
omitted. The 23 sub determinant is generically of order 7]. The mass matrix in (jUjl 
does not describe the most general perturbation of the zeroth order texture (|40|) . We 
have implicitly assumed a symmetry between and u T which is preserved by the 
perturbations, at least at the level of the order of magnitudes. The perturbed texture 
can also arise when the zeroes of the lopsided Dirac matrix in (|3~5j) are replaced 



by small quantities. It is possible to construct models along this line based on a 
spontaneously broken U(l)p flavor symmetry, where 5, e and 77 are given by positive 
powers of one or more symmetry breaking parameters. Moreover, by playing with 
the U(1 )f charges, we can adjust, to certain extent, the relative hierarchy between 77, 
e and 5 60 61 1 63 1 64 1 65 1 66 ^ ag we w jjj gee j n sec ^i on g_ The texture (pHj) can also be 

generated in SUSY models with i?-parity violation!^. 

Let us come back to the mass matrix m v of eq. ()44j1 . After a first rotation by an 
angle #23 close to 7r/4 and a second rotation with 6*13 ~ e, we get 

/ 5 + e 2 e \ 
m v Km e ^0 , (45) 

V 02/ 

up to order one coefficients in the small entries. To obtain a large solar mixing angle, 
we need \q — 8\ < e. In realistic models there is no reason for a cancellation between 
independent perturbations and thus we assume both 5 < e and 77 < e. 

Consider first the case 5 ~ e and 77 < e. The solar mixing angle 9 12 is large 
but not maximal, as indicated by the LA solution. We also have Am 2 atm w 4m 2 , 



Am Ln » Am lt m t 2 and 



m ee ^ Am 2 sun . (46) 

If 77 ~ e and 5 e, we still have a large solar mixing angle and Am 2 un ~ e 2 Am 2 . tm , 
as before. However m PP will be much smaller than the estimate in (1461). This is the 



case of the models based on the above mentioned U(l)p flavor symmetry that, at 
least in its simplest realization, tends to predict 5 ~ e 2 . In this class of models we 
find 



/ Am 



2 



m ee « ^Am 2 un j , (47) 

below the sensitivity of the next generation of planned experiments. It is worth to 
mention that in both cases discussed above, we have 



which might be close to the present experimental limit if the oscillation frequency of 
the LA solution for solar neutrinos is in the upper part of its allowed range. 

If both 5 and 77 are much smaller than e, the 12 block of m u has an approximate 
pseudo-Dirac structure and the angle #12 becomes maximal. This situation is typical 
of some models where leptons have U(l)p charges of both signs whereas the order 
parameters of U(l)p breaking have all charges of the same sign^^. We have two 
eigenvalues approximately given by ±m e. As an example, we consider the case 
where 77 = and 5 « e 2 . We find sin 2 29 u ~ 1 — e 2 /4, Am 2 sun ps m 2 e 3 



In order to recover the LA solution we would need a relatively large value of e. 
But this is in general problematic because, on the one hand the presence of a large 
perturbation raises doubts about the consistency of the whole approach and, on the 
other hand, in existing models where all fermion sectors are related to each other, e 
is never larger than the Cabibbo angle. 

Summarising, within normal hierachical models there is enough flexibility to re- 
produce in a natural way the experimental frequencies and mixing angles. In partic- 
ular the lopsided matrix solution of the large atmospheric mixing, inspired by SU(5), 
where the large atmospheric mixing arises from the charged lepton sector, can be 
extended rather naturally to also account for the solar sector and for the small #13 
mixing angle. 



8.1. Semi- anarchy 

We have seen that anarchy is the absence of structure in the neutrino sector. Here 
we consider an attenuation of anarchy where the absence of structure is limited to 
the 23 sector. The typical texture is in this case: 

/ 5 e e \ 

m u «m ell , (50) 
\e 1 l) 

where 5 and e are small and by 1 we mean entries of 0(1) and also the 23 determinant 
is of 0(1). We see that this texture is similar to eq. (f4*4*|) when r) ~ 0(1). Clearly, 
in general we would expect two mass eigenvalues of order 1, in units of m, and one 
small, of order 5 or e 2 . This pattern does not fit the observed solar and atmospheric 
observed frequencies. However, given that the ratio r = (Am^J^/Am^ ~ 1/40 
is not too small, we can assume that the small value of r is generated accidentally, 
as for anarchy. We see that, if we proceed with the same change of basis as from 
eq. ()44|) to eq. ()45|) . it is sufficient that by chance 77 ~ 5 + e 2 in order to obtain the 
correct value of r with large #23 and 6*12 and small #13 ~ e. The natural smallness of 
#13 is the main advantage over anarchy. We will come back to this class of models in 
a following section. 



9. Grand Unified Models of Fermion Masses 



We have seen that the smallness of neutrino masses interpreted via the see-saw 
mechanism directly leads to a scale A for L non- conservation which is remarkably 
close to Mqut- Thus neutrino masses and mixings should find a natural context 
in a GUT treatment of all fermion masses. The hierarchical pattern of quark and 
lepton masses, within a generation and across generations, requires some dynamical 
suppression mechanism that acts differently on the various particles. This hierarchy 
can be generated by a number of operators of different dimensions suppressed by 



inverse powers of the cut-off A c of the theory. In some realizations, the different powers 
of 1/A C correspond to different orders in some symmetry breaking parameter Vf arising 
from the spontaneous breaking of a flavour symmetry. In the next subsections we 
describe some simplest models based on SU(5) x U(1)f and on SO(10) which illustrate 
these possibilities'^. It is notoriously difficult to turn these models into fully realistic 
theories, due to well-known problems such as the doublet-triplet splitting, the proton 
lifetime, the gauge coupling unification beyond leading order and the wrong mass 
relations for charged fermions of the first two generations. Some of these problems 
can be solved by adopting the elegant idea of GUT's in extra dimensions"^. Here 
we adopt the GUT framework simply as a convenient testing ground for different 
neutrino mass scenarios. 

9.1. Models Based on Horizontal Abelian Charges 

We discuss here some explicit examples of grand unified models in the framework 
of a unified SUSY SU(5) theory with an additional U(l)p flavour symmetry. The 
SU(5) generators act "vertically" inside one generation, while the U(1)f charges are 
different "horizontally" from one generation to the other. If, for a given interaction 
vertex, the U(l)p charges do not add to zero, the vertex is forbidden in the symmetric 
limit. But the symmetry is spontaneously broken by the VEV's Vf of a number of 
"flavon" fields with non-vanishing charge. Then a forbidden coupling is rescued but 
is suppressed by powers of the small parameters Vf/A c with the exponents larger for 
larger charge mismatch^!. We expect M GUT Vf ^ A c ^ M P \. Here we discuss 
some aspects of the description of fermion masses in this framework. 

In these models the known generations of quarks and leptons are contained in 
triplets and (i = 1,2,3) corresponding to the 3 generations, transforming 
as 10 and 5 of SU(5), respectively. Three more SU(5) singlets ty\ describe the RH 
neutrinos. In SUSY models we have two Higgs multiplets H u and H^, which transform 
as 5 and 5 in the minimal model. The two Higgs multiplets may have the same or 
different charges. In all the models that we discuss the large atmospheric mixing 
angle is described by assigning equal flavour charge to muon and tau neutrinos and 
their weak SU(2) partners (all belonging to the 5 = (l,d c ) representation of SU(5)). 
Instead, the solar neutrino oscillations can be obtained with different, inequivalent 
charge assignments. There are many variants of these models: fermion charges can all 
be non-negative with only negatively charged flavons, or there can be fermion charges 
of different signs with either flavons of both charges or only flavons of one charge. 
We can have that only the top quark mass is allowed in the symmetric limit, or that 
also other third generation fermion masses are allowed. The Higgs charges can be 
equal, in particular both vanishing or can be different. We can arrange that all the 



structure is in charged fermion masses while neutrinos are anarchical. 



9.1.1. F(fermions)> 

Consider, for example, a simple model with all charges of matter fields being non- 
negative and containing one single flavon 9 of charge F= —1. For a maximum of 
simplicity we also assume that all the third generation masses are directly allowed 
in the symmetric limit. This is realized by taking vanishing charges for the Higgses 
and for the third generation components ^3°, ^3 and ^>\. If we define F(\l/f ) = qf- 
(R = 10, 5, 1; i = 1, 2, 3), then the generic mass matrix m has the form 



m 



(yu\ q ? +q ? yi2\ q ? +q °' yi3\ q ? +q ?'\ 

V2l\ q * +q ? V22\ q * +q *' y23\ q ? +q ?' 

\y 31 \^+«?' yaa \<#+£ y^\ q * +q *' ) 



(51) 



where all the i/ij are dimensionless complex coefficients of order one and m u , = mf, 
m D and M arise by choosing (R, R') = (10, 10), (5, 10), (1, 5) and (1, 1), respectively. 
We have A = (0) / A c and the quantity v represents the appropriate VEV or mass 
parameter. The models with all non-negative charges and one single flavon have 
particularly simple factorization properties. For instance in the see-saw expression for 
w„ = rn^M^mo the dependence on the qj charges drops out and only that from qf 
remains. In addition, for the neutrino mixing matrix C/^, which is determined by m v 
in the basis where the charged leptons are diagonal, one can prove that Uij ~ A' 9i_ ^ ', 
in terms of the differences of the 5 charges, when terms that are down by powers 
of the small parameter A are neglected. Similarly the CKM matrix elements are 
approximately determined by only the 10 charges V£j KM ~ A' 9 * "^ '. If the 
symmetry breaking parameter A is numerically close to the Cabibbo angle, we can 
choose: 

(ql°,q 1 2°,ql ) = (3,2,0) , (52) 

thus reproducing V us ~ A, V c b ~ A 2 and V u b ~ A 3 . The same q}° charges also fix 
m u : m c : m t ~ A 6 : A 4 : 1. The experimental value of m u (the relevant mass values 
are those at the GUT scale: m = m^Moir) would rather prefer q\° = 4. Taking 
into account this indication and the presence of the unknown coefficients ~ 0(1) 
it is difficult to decide between q\° = 3 or 4 and both are acceptable. Of course the 
charges (q\°, ql°, ql°) = (2, 1, 0) would represent an equally good choice, provided we 
appropriately rescale the expansion parameter A. Turning to the 5 charges, if we take 

(qlqlq!) = (b,o,o) b>o , (53) 

together with eq. (|52"|) we get the patterns : m s : m;, ~ m e : m M : m T ~ X 3+b : A 2 : 
1. Moreover, the 22, 23, 32, 33 entries of the effective light neutrino mass matrix m u 
are all 0(1), thus accommodating the nearly maximal value of S23. The small non 



diagonal terms of the charged lepton mass matrix cannot change this. We obtain, 
where arbitrary 0(1) coefficients are omitted: 




m v = \ b 1 lU (A,SA) , (54) 



where v u is the VEVs of the Higgs doublet giving mass to the up quarks and all the 
entries are specified up to order one coefficients. If we take v u ~ 250 GeV, the mass 
scale A of the heavy Majorana neutrinos turns out to be close to the unification scale, 
A ~ 10 15 GeV. 

If b vanishes, then the light neutrino mass matrix will be structure-less and we 
recover the anarchical (A) picture of neutrinos discussed in section 6.1. In a large 
sample of anarchical models, generated with random coefficients, the resulting neu- 
trino mass spectrum can exhibit either normal or inverse hierarchy. For down quarks 
and charged leptons we obtain a weakened hierarchy, essentially the square root than 
that of up quarks. 

If b is positive, then the light neutrino mass matrix will be structure-less only in 
the (2,3) sub-sector and we get semi-anarchical (SA) models, introduced in section 
7.2. In this case, the neutrino mass spectrum has normal hierarchy. However, unless 
the (2,3) sub-determinant is accidentally suppressed, atmospheric and solar oscillation 
frequencies are expected to be of the same order and, in addition, the preferred solar 
mixing angle is small. Nevertheless, such a suppression can occur in a fraction of semi- 
anarchical models generated with random, order one coefficients. The real advantage 
over the fully anarchical scheme is represented by the suppression in U e ^. 

Note that in all previous cases we could add a constant to qf, for example by 
taking (qf, q\, q\ ) = (2 + 6,2,2). This would only have the consequence to leave 
the top quark as the only unsuppressed mass and to decrease the resulting value of 
tan (3 = v u /vd down to X 2 mt/mi ) . A constant shift of the charges q\ might also provide 
a suppression of the leading v c mass eigenvalue, from A c down to the appropriate scale 
A. One can also consider models where the 5 and 5 Higgs charges are different, as in 
the "realistic" SU(5) model of ref. Also in these models the top mass could be 
the only one to be non- vanishing in the symmetric limit and the value of tan j3 can 
be adjusted. 



9.1.2. F(fermions) and F(flavons) of both signs 

Models with naturally large 23 splittings are obtained if we allow negative charges 
and, at the same time, either introduce flavons of opposite charges or stipulate that 
matrix elements with overall negative charge are put to zero. For example, we can 
assign to the fermion fields the set of F charges given by: 



(ql ^?,^ ) = (3,2,0) 



(g?,<idb = (6,0,0) 6>2a>0 



(55) 



We consider the Yukawa coupling allowed by U(l)F-neutral Higgs multiplets in the 5 
and 5 SU(5) representations and by a pair 9 and 9 of SU(5) singlets with F= 1 and 
F= —1, respectively. If b = 2 or 3, the up, down and charged lepton sectors are not 
essentially different than in the SA case. Also in this case the 0(1) off-diagonal entry 
of mi, typical of lopsided models, gives rise to a large LH mixing in the 23 block which 
corresponds to a large RH mixing in the d mass matrix. In the neutrino sector, after 
diagonalization of the charged lepton sector and after integrating out the heavy RH 
neutrinos we obtain the following neutrino mass matrix in the low-energy effective 
theory: 

f\ 2b X b X b \ 2 

m v = X b l + X a X' a l + X a X' a ^ (H), (56) 
V X b l + A a A /a l + A a A' a / 

where A' is given by (9) / A c and A as before denotes the large mass scale associated 
to the RH neutrinos: A ^> v U)d . The 0(1) elements in the 23 block are produced 
by combining the large LH mixing induced by the charged lepton sector and the 
large LH mixing in mo- A crucial property of m v is that, as a result of the see-saw 
mechanism and of the specific U(1)f charge assignment, the determinant of the 23 
block is automatically of 0(A a A /a ) (for this the presenc e of negative charge values, 
leading to the presence of both A and A' is essential EaESu ). The neutrino mass matrix 
of eq. (J5T)|) is a particular case of the more general pattern presented in eq. (I44JI . for 
5 « X 2b , f k ^ and rj « A a A' a . If we take A « A', it is easy to verify that the 
eigenvalues of m u satisfy the relations: 



mi : m2 '■ 1^3 



X 2(b-a) . X 2a. 1 _ ^ 



The atmospheric neutrino oscillations require m 2 ~ 10 -3 eV 2 . The squared mass 
difference between the lightest states is of 0(X 4a ) m|, not far from the LA solution 
to the solar neutrino problem if we choose a = 1. In general U e ^ is non- vanishing, 
of 0(X b ). Finally, beyond the large mixing in the 23 sector, m v provides a mixing 
angle Q\i ~ \ b - 2a j n ^he 12 sector. When b = 2a, as for instance in the case 6 = 2 
and a — 1, the LA solution can be reproduced and the resulting neutrino spectrum 
is hierarchical (H). 

Alternatively, an inversely hierarchical (IH) spectrum can be obtained by choosing: 

(Si 10 ,^ ) = (3,2,0) 

{q\,ql,ql) = (1,-1,-1) 

(?i, 92,53) = (-1,1,0) 

(<Hf u ,q Hd ) = (0,1) . (58) 



Due to the non- vanishing charge of the Ha Higgs doublet, in the charged lepton sector 
we recover the same pattern previously discussed. The light neutrino mass matrix is 
given by: 

/ A 2 1 1 
m v = 1 A' 2 A' 2 
V 1 A' 2 A' 2 

The ratio between the solar and atmospheric oscillation frequencies is not directly 
related to the sub-determinant of the block 23, in this case. 



(IH) 



(59) 



Table 2: Models and their flavour charges. 



Model 


*10 


*5 




(H u , H d ) 


Anarchical (A) 


(3,2,0) 


(0,0,0) 


(0,0,0) 


(0,0) 


Semi-Anarchical (SA) 


(2,1,0) 


(1,0,0) 


(2,1,0) 


(0,0) 


Hierarchical (Hi) 


(6,4,0) 


(2,0,0) 


(1,-1,0) 


(0,0) 


Hierarchical (Ha) 


(5,3,0) 


(2,0,0) 


(1,-1,0) 


(0,0) 


Inversely Hierarchical (IHj) 


(3,2,0) 


(1,-1,-1) 


(-i,+i,o) 


(o,+i) 


Inversely Hierarchical (IHjj) 


(6,4,0) 


(1,-1,-1) 


(-i,+i,o) 


(o,+i) 



A representative set of models is listed in table 2. Note that in some cases the 
charges for ty w have been changed from (3,2,0) (our reference values in eqs. (|52|) . 
(JHSJ), and (jSHI)) to (6,4,0) or (5,3,0). These values are a posteriori better suited to 
reproduce the moderate level of hierarchy implied by the present neutrino oscillation 
data. Since the neutrino mixing parameters are completely independent on the 10 
charges, this change is only important for a better fit to quark and charged lepton 
masses and mixings once a rather large value of A is derived from the neutrino data. 
The hierarchical and the inversely hierarchical models may come into several varieties 



depending on the number and the charge of the flavour symmetry breaking (FSB) 
parameters. Above we have considered the case of two (II) oppositely charged flavons 
with symmetry breaking parameters A and A'. It may be noticed that the presence of 
two multiplets 9 and 9 with opposite F charges could hardly be reconciled, without 
adding extra structure to the model, with a large common VEV for these fields, due 
to possible analytic terms of the kind (99) n in the superpotential. Therefore it is 
instructive to explore the consequences of allowing only the negatively charged 9 field 
in the theory, case I. In case I, it is impossible to compensate negative F charges in 
the Yukawa couplings and the corresponding entries in the neutrino mass matrices 
vanish. Eventually these zeroes are filled by small contributions, arising, for instance, 
from the diagonalization of the charged lepton sector or from the transformations 
needed to make the kinetic terms canonical. 

Another important ingredient is represented by the see-saw mechanism^. Hier- 
archical models and semi- anarchical models have similar charges in the (10, 5) sectors 
and, in the absence of the see-saw mechanism, they would give rise to similar results. 
Even when the results are expected to be independent from the charges of the RH 
neutrinos, as it is the case for the anarchical and semi-anarchical models, the see-saw 
mechanism can induce some sizeable effect in a statistical analysis. For this reason, 
for each type of model, but the normal-hierarchical ones (the mechanism for the 23 
sub-determinant suppression is in fact based on the see-saw mechanism), it is inter- 
esting to study the case where RH neutrinos are present and the see-saw contribution 
is the dominant one (SS) and the case where they are absent and the mass matrix is 
saturated by the non-renormalizable contribution (NOSS). 

With this classification in mind, we can distinguish the following type of models, 
all supported by specific choices of U(l) charges: Ass, Anoss, SAss, SAnoss, H(ss,i), 

H(ss,ii)> IH(ss,i), IH(ss,ii), IH (NO ss,i) and IH( NO ss,ii)- 

It is interesting to quantify the ability of each model in reproducing the observed 
oscillation parameters. For anarchy, it has been observed that random generated, 
order-one entries of the neutrino mass matrices (in appropriate units), correctly fit 
the experimental data with a success rate of few percent. It is natural to extend this 
analysis to include also the other models based on SU(5) x U(l)^^, which have mass 
matrix elements defined up to order-one dimensionless coefficients y$j (see eq. IS*!]) . For 
each model, successful points in parameter space are selected by asking that the four 
observable quantities 0\ = r = Am 2 2 /| Am 2 3 |, 0% = tan 2 9u, O3 = \U e s\ = | sin#i 3 | 



and O4 = tan #23 fall in the approximately 3a allowed ranges 1 1 s 

0.018 < r < 0.053 
\U e3 \ < 0.23 
0.30 < tan 2 12 < 0.64 
0.45 < tan 2 23 < 2.57 



2n (6°) 



The coefficients of the neutrino sector are random complex numbers with absolute 
values and phases uniformly distributed in intervals I = [0.5, 2] and [0, 27r] respec- 



tively. The dependence of the results on these choices can be estimated by varying X. 
For each model an optimization procedure selects the value of the flavour symmetry 
breaking parameter A = A' that maximizes the success rate. The success rates are dis- 
played in figs. El and HI separately for the SS and NOSS cases. The two sets of models 
have been individually normalized to give a total rate 100. From the histograms in 
figs. ^andEJwe see that normal hierarchy models are neatly preferred over anarchy 
and inverse hierarchy in the context of these SU(5)xU(l) models. In particular, in 
the SS case, the Hit models with normal hierarchy, two oppositely charged flavons 
and suppressed 23 sub-determinant are clearly preferred. Models of the type Hi are 
disfavoured. For the relatively large values of the expansion parameter required to fit 
r, they tend to predict too large \U e s\ and tan 2 9 12 > 1. We recall that for the chosen 
charge values the H n model is of the lopsided type. In the NOSS case the see-saw 
suppression of the 23 determinant is clearly not operative and all normal hierarchy 
models coincide with SA. 

70 
60 
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30 
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Figure 1: Relative success rates for the LA solution, with see-saw. The sum of the rates has been 
normalized to 100. The results correspond to the default choice X = [0.5,2], and to the following 
values of A = A': 0.2, 0.25, 0.35, 0.45, 0.45, 0.25 for the models A ss , SA SS , H (SS j I) , H( SS)I) , IH(ss,n) 
and 111(33,1), respectively. The error bars represent the linear sum of the systematic error due to the 
choice of X and the statistical error (see text). 

An interesting question is whether the disfavouring of IH models that we find 
in our SU(5)xU(l) framework can be extended to a more general context. In the 
limit of vanishing A and A' the IH texture (see eq. (J59)l) becomes close to that of 
bimaximal mixing and #13 = (actually with r = 0). In our U(l) models r rs [L^l ~ 
I tan 2 6*12 — 1 1 ~ 0(A 2 ) (for A = A'). In particular the charged lepton mixings cannot 
displace too much 9yi from its maximal value because the small value of the electron 
mass forces a sufficiently large value of the relevant charges, which in turn implies 
that the charged lepton mixing correction to 9\2 is small. We have already mentioned 
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Figure 2: Relative success rates for the LA solution, without see-saw. The sum of the rates has been 
normalized to 100. The results correspond to the default choice T = [0.5,2], and to the following 
values of A = A': 0.2, 0.2, 0.25, 0.25 for the models A NO ss, SA NO ss, IH(nossji), and IH (NO ss,i), 
respectively (in our notation there are no H(noss,i)> H(nossji) models). The error bars represent 
the linear sum of the systematic error due to the choice of X and the statistical error (see text) . 

that corrections from the charged lepton sector can in principle bring the predictions 
of a neutrino matrix of the bimixing type in agreement with the data and that the 
smallness of S13 induces strong constraints. In the particular setup of U(l)p models 
we have seen that charged lepton corrections are too small to make the solar angle 
sufficiently different from maximal. 

Leptonic mixing in SUSY GUTs with RH neutrinos are potential sources of lepton 
flavour violating (LFV) processes beyond neutrino oscillations^^. The observable 
effects are difficult to estimate since they are sensitive to both the details of the 
SUSY breaking mechanism and to the low-energy superparticle spectrum. In models 
with gravity mediated SUSY breaking and universal boundary conditions for the soft 
breaking terms at the cut-off scale A, running effects give rise to off-diagonal slepton 
masses Smfj LL proportional, at leading order, to = (yl\og(A/M)y u )ij. Current 
bounds on LFV transitions U — > Iff translate into an upper bound on the combination 
\Cij\, depending on tan/? and soft mass parameters. If, for instance, tan/3 = 10, the 
present experimental bound on BR(fj, — > ej) already excludes C 2 i > 1 for the most 
plausible values of slepton and gaugino masses. In the models considered in this 
section C21 is dominated by the couplings to the RH neutrino of the third generation 
and, assuming (A/M 3 ) « O(100), we roughtly expect 1 2/322/31 1 < 0-2- This constraint 
is most easily respected by models with inverse hierarchy. For semianarchy and 
normal hierarchy the constraint is almost saturated, whereas anarchy tends to violate 
it. Future improvements in the experimental sensitivity could lead to a significant 
selection of the models. 



In conclusion, models based on SU(5) x U(l)p are clearly toy models that can 
only aim at a semiquantitative description of fermion masses. In fact only the order 
of magnitude of each matrix entry can be specified. However it is rather impressive 
that a reasonable description of fermion masses, now also including neutrino masses 
and mixings, can be obtained in this simple context, which is suggestive of a deeper 
relation between gauge and flavour quantum numbers. There are 12 mass eigenvalues 
and 6 mixing angles that are specified, modulo coefficients of order 1, in terms of 
a bunch of integer numbers (from half a dozen to a dozen), the charges, plus 1 
or more scale parameters. Moreover all possible type of mass hierarchies can be 
reproduced within this framework. In a statistically based comparison, the range 
of r and the small upper limit on U e s are sufficiently constraining to make anarchy 
neatly disfavoured with respect to models with built-in hierarchy. If only neutrinos 
are considered, one might counterargue that hierarchical models have at least one 
more parameter than anarchy, in our case the parameter A. However, if one looks 
at quarks and leptons together, as in the GUT models that we consider, then the 
same parameter that plays the role of an order parameter for the CKM matrix, for 
example, the Cabibbo angle, can be successfully used to reproduce also the hierarchy 
implied by the present neutrino data. 

9.2. GUT Models based on SO(10) 

The fermion sector of SO (10) grand unified theories has remarkable properties. It 
is automatically anomaly-free, independently from the representation content. The 
RH neutrino provides the completion of a SM family into a 16 representation, thus 
offering a natural ground for the see-saw mechanism and baryogenesis through lep- 
togenesis. Moreover, the scale of lepton number violation, determined by the gauge 
symmetry breaking pattern, can be smaller than the cut-off A c of the theory, and 
closer to the grand unification scale itself as suggested by the mass scale associ- 
ated to atmospheric neutrino oscillations. 

In their simplest realizations SO(10) models are left-right symmetric at the GUT 
scale and we would expect similar mixing angles for quarks and leptons of both 
chiralities. In left-right symmetric models smallness of left mixings implies that also 
right-handed mixings are small, so that all mixings tend to be small, unless non- 
renormalizable mass operators with a suitable flavour pattern are introduced. In 
such a context, accommodating the observed mixing properties of quarks and leptons 
appears more problematic than in SU(5) theories, at first sight. 

One possibility is to exploit the see-saw mechanism to enhance the light neutrino 
mixing angles. As illustrated in eq. (}38j) in a 2 by 2 context, in order to have large 
or even maximal mixing in m u , we do not necessarily need large mixing angles in 
and M. We can obtain a large mixing starting from nearly diagonal mass matrices 
m D and M, provided a RH neutrino, equally coupled to and z/ r , is sufficiently light 



An other possibility is to argue that perhaps what appears to be large is not 
that large after all. The typical small parameter that appears in the mass matrices 
is A ~ \Jmd/m s ~ y^vT^V ~ 0-20 — 0.25. This small parameter is not so small 
that it cannot become large due to some peculiar accidental enhancement: either a 
coefficient of order 3, or an exponent of the mass ratio which is less than 1/2 (due for 
example to a suitable charge assignment), or the addition in phase of an angle from 
the diagonalization of charged leptons and an angle from neutrino mixing. Typically, 
by exploiting the freedom in the parameter space, in this set of models a large #23 
may be accommodated. The large mixing angle for solar neutrinos requires how- 
ever the introduction of ad hoc terms, like for instance higher-dimensional operators 

contribu ting to the light neutrino masses independently from the see-saw mechanism 
84 85 86, 

Alternatively, to avoid the introduction of ad hoc non-renormalizable operators, it 
is possible to enlarge the Higgs content and consider, for example, an SO (10) model 
where all fermion mass matrices originate from renormalizable interactions of matter 
fields in three 16 representations with two Higgs multiplets, a 10// and a 126//^^: 

C Y = 10h 16 2/10I6 + 126// 16 y 126 16 , (61) 

where yio and yi 2 e are two symmetric matrices in flavour space. Both yio and yi 2 Q 
contribute to the Dirac mass matrices, with the characteristic factor -3 for the yi 2 o 
between the (u, d) and (u, e) sectors: 

m d = a yio + (3 yi 26 

m e = a y w - 3p y 126 , (62) 

and the correct mass relations for first and second generations can be accommodated. 
In the most general case, when 126// acquires VEVs in its SU(2)i singlet and triplet 
components, both RH and LH Majorana masses can arise in the neutrino sector. 
They are both proportional to yi 26 and give rise to type I and type II see-saw, re- 
spectively. By assuming dominance of the type II contribution, we find an interesting 
link between large atmospheric mixing angle and b — r unification . In this case 
the light neutrino mass matrix is proportional to 1/126 and, from eq. ()62|). it can be 
directly related to m e and m d : 

m v oc ma — m e . (63) 
If both m e and md have the approximate pattern 

'l *) (64) 



in the 23 sector, then b — r unification forces a cancellation in the 33 entry of m v , 
thus allowing for a large 23 neutrino mixing angle. Most of the Yukawa parameters 



can be determined by the quark masses and mixing angles and from the charged 
lepton masses. This allows to predict the atmospheric and solar mixing parameters 
within a range which is still experimentally allowed and \U e s\ ~ 0.16, not far from the 
present upper bound A drawback of the model is the occurrence of two pairs of 
Higgs doublets with non-vanishing VEVs, of which only one combination is allowed 
to remain light. This clearly makes the doublet-triplet splitting problem even more 
complicated than in minimal models. 

Finally, we can abandon the idea that the model is left-right symmetric at the GUT 
scale. In this case, the mechanism discussed in section 8, based on asymmetric mass 
mat rices, can be em bedded in an SO(10) grand-unified theory in a rather economic 
way IBISZISIEDISII The 33 entries of the fermion mass matrices can be obtained 
through the coupling 16316310^ among the fermions in the third generation, I63, 
and a Higgs tenplet 10h- The two independent VEVs of the tenplet v u and Vd give 
mass, respectively, tot/v T and b/r. The key point to obtain an asymmetric texture is 
the introduction of an operator of the kind 16216^16316^ . This operator is thought 
to arise by integrating out an heavy 10 that couples both to I62I6// and to 16316^. 
If the IQh develops a VEV breaking SO(10) down to SU(5) at a large scale, then, 
in terms of SU(5) representations, we get an effective coupling of the kind 52IO3BH, 
with a coefficient that can be of order one. This coupling contributes to the 23 entry 
of the down quark mass matrix and to the 32 entry of the charged lepton mass matrix, 
realizing the desired asymmetry. To distinguish the lepton and quark sectors one can 
further introduce an operator of the form 16il6jl0#45#, = 2, 3), with the VEV 
of the 45 H pointing in the B — L direction. Additional operators, still of the type 
16jl6jl6//16^ can contribute to the matrix elements of the first generation. The 
mass matrices look like: 
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where r] 8, 8' e 1 and a, b and a are of order 0(1). In the charged fermion 
sector, the parameter rj, 8, 8', e and a are determined from the lepton masses and 
from a subset of quark masses and mixing angles, leading to six successful predictions. 
In the neutrino sector, the lopsidedness of mi is responsible for the large atmospheric 
mixing angle, whereas the parameters a and b can be adjusted to obtain the solar 
mixing angle and the ratio between solar and atmospheric squared mass differences. 



The model predicts small values for \U e s\, in a range accessible only to future neutrino 
factories ^ . 

Models based on SO(10) times a flavour symmetry are more difficult to construct 
because a whole generation is contained in the 16, so that, for example for U(l)p, one 
would have the same value of the charge for all quarks and leptons of each generation, 
which is too rigid. This problem can be circumvented if not all the observed fermions 
in a given generation belong to a single 16 representation!^. 



10. Conclusion 



By now there are rather convincing experimental indications for neutrino oscil- 
lations. The direct implication of these findings is that neutrino masses are not all 
vanishing. As a consequence, the phenomenology of neutrino masses and mixings is 
brought to the forefront. This is a very interesting subject in many respects. It is 
a window on the physics of GUTs in that the extreme smallness of neutrino masses 
can only be explained in a natural way if lepton number conservation is violated. If 
so, neutrino masses are inversely proportional to the large scale where lepton number 
is violated. Also, the pattern of neutrino masses and mixings interpreted in a GUT 
framework can provide new clues on the long standing problem of understanding the 
origin of the hierarchical structure of quark and lepton mass matrices. 

Neutrino oscillations only determine differences of mf values and the actual scale 
of neutrino masses remain to be experimentally fixed. The detection of Ovflfl decay 
would be extremely important for the determination of the overall scale of neutrino 
masses, the confirmation of their Majorana nature and the experimental clarification 
of the ordering of levels in the associated spectrum. The recent results from cosmol- 
ogy indicate that neutrino masses are not a major fraction of the cosmological mass 
density Q u ^ 1.5%. The decay of heavy right-handed neutrinos with lepton number 
non-conservation can provide a viable and attractive model of baryogenesis through 
leptogenesis. The measured oscillation frequencies and mixings are remarkably con- 
sistent with this attractive possibility. 

While the existence of oscillations appears to be on a ground of increasing solidity, 
many important experimental challenges remain. For atmospheric neutrino oscilla- 
tions the completion of the K2K experiment, which was delayed by the accident that 
has seriously damaged the Superkamiokande detector, is important for a terrestrial 
confirmation of the effect and for an independent measurement of the associated pa- 
rameters. In the near future the experimental study of atmospheric neutrinos will be 
further pursued with long baseline measurements by MINOS, OPERA, ICARUS. For 
solar neutrinos the continuation of SNO, KamLAND and the data from Borexino will 
lead to a more precise determination of the parameters of the LA solution. Finally a 
clarification by MINIBOONE of the issue of the LSND alleged signal is necessary, in 
order to know if 3 light neutrinos are sufficient or additional sterile neutrinos must be 



introduced, in spite of the apparent lack of independent evidence in the data for such 
sterile neutrinos and of the fact that attempts of constructing plausible and natural 
theoretical models have not led so far to compelling results. Further in the future 
there are projects for neutrino factories and/or superbeams aimed at precision mea- 
surements of the oscillation parameters and possibly the detection of CP violation 
effects in the neutrino sector. 

Pending the solution of the existing experimental ambiguities a variety of theoret- 
ical models of neutrino masses and mixings are still conceivable. Among 3-neutrino 
models we have described a number of possibilities based on degenerate, inverted hi- 
erarchy and normal hierarchy type of spectra. The normal hierarchy option appears 
to us as the most straightforward and flexible framework. In particular the large 
atmospheric mixing can arise from lopsided matrices. Then the observed frequencies 
and the large solar angle can also be obtained without fine tuning in models where 
the 23 subdeterminant is automatically suppressed. 

The fact that some neutrino mixing angles are large and even nearly maximal, 
while surprising at the start, was eventually found to be well compatible with a 
unified picture of quark and lepton masses within GUTs. The symmetry group at 
Mqut could be either (SUSY) SU(5) or SO(10) or a larger group. For example, we 
have seen that models based on anarchy, semianarchy, inverted hierarchy or normal 
hierarchy can all be naturally implemented by simple assignments of U(l)p horizontal 
charges in a semiquantitative unified description of all quark and lepton masses in 
SUSY SU(5)x U(1)f- Actually, in this context, if one adopts a statistical criterium, 
hierarchical models appear to be preferred over anarchy and among them normal 
hierarchy appears the most likely. Note that in almost all of the existing models of 
neutrino mixings the atmospheric angle is large but not maximal. If it experimentally 
turned out that indeed #23 i s maxim al with good accuracy then very special classes 
of models would be selected fc6 | 43 | 94) 

All we know about neutrino masses is well in harmony with the idea and the mass 
scale of GUT's. As a consequence neutrino masses have added phenomenological 
support to this beautiful idea and to the models of physics beyond the Standard 
Model that are compatible with it. In particular, if we consider the main classes of 
new physics that are currently contemplated, like supersymmetry, technicolour, large 
extra dimensions at the TeV scale, little Higgs models etc, it is clear that the first one 
is the most directly related to GUT's. SUSY offers a well defined model computable 
up to the GUT scale and is actually supported by the quantitative success of coupling 
unification in SUSY GUT's. For the other examples quoted all contact with GUT's 
is lost or at least is much more remote. In this sense neutrino masses fit particularly 
well in the SUSY picture that sofar remains the standard way beyond the Standard 
Model. 
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